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We study a frictionless pendulum subject to multiplicative random noise. Because of destructive 
interference between the angular displacement of the system and the noise term, the energy fluctua- 
tions are reduced when the noise has a non-zero correlation time. We derive the long time behavior 
of the pendulum in the case of Ornstein-Uhlenbeck noise by a recursive adiabatic elimination pro- 
cedure. An analytical expression for the asymptotic probability distribution function of the energy 
is obtained and the results agree with numerical simulations. Lastly, we compare our method to 
other approximation schemes. 
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I. INTRODUCTION 



The behavior of a nonlinear dynamical system is strongly modified when randomness is taken into account: noise 
can shift bifurcation thresholds, create new phases (noise-induced transitions), or even generate spatial patterns 
US Hll 13 The interplay of noise with nonlinearity gives rise to a variety of phenomena that constantly motivate 
new research of theoretical and practical significance. Stochastic resonance [3 and biomolecular Brownian motors 
0,0 are celebrated examples of nonlinear random systems of current interest. In particular, stochastic ratchets have 
generated a renewed interest in the study of simple mechanical systems subject to random interactions, the common 
ancestor of such models being Langevin's description of Brownian motion. Many unexpected phenomena appear when 
one generalizes Langevin's equations to include, e.g., inertial terms, nonlinearities, external (multiplicative) noise or 
noise with finite correlation time; each of these new features opens a field of investigations that calls for specific 
techniques or approximation schemes _;9|. Nonlinear oscillators with parametric noise are often used as paradigms 
for the study of these various effects and related mathematical methods ^} . The advantage of such models is that 
they have an appealing physical interpretation and appear as building blocks in many different fields; they can be 
simulated on a computer or constructed as real electronic or mechanical systems Il2j . Moreover, the mathematical 
apparatus needed to analyze them remains relatively elementary (as compared to the perturbative field-theoretical 
methods required for spatio-temporal systems and can be expected to yield exact and rigorous results. 

In the present work, we study the motion of a frictionless pendulum with parametric noise, which can be physically 
interpreted as a randomly vibrating suspension axis. We show that the long time behavior of a stochastic pendulum 
driven by a colored noise with finite correlation time is drastically different from that of a pendulum subject to white 
noise. Whereas the average energy of the white-noise pendulum is a linear function of time, that of the colored-noise 
pendulum grows only as the the square-root of time. Our analysis is based on a generalization of the averaging 
technique that we have used previouly for nonlinear oscillators subject to white noise: in [bH IT5I fill an effective 
dynamics for the action variable of the system is derived after integrating out the fast angular variable, and is then 
exactly solved. However, this averaging technique as such ceases to apply to systems with colored noise when the 
time scale of the fast variable becomes smaller than the correlation time of the noise. Correlations between the 
fast variable and the noise modify the long time scaling behavior of the system and therefore must be taken into 
account. We shall develop here a method that systematically retains these correlation terms before the fast variable is 
averaged out. This will allow us to derive analytical expressions for the asymptotic probability distribution function 
(P.D.F.) of the energy of the stochastic pendulum, and to deduce the long time behavior of the system. Our analytical 
results are verified by numerical simulations. Finally, we shall compare our method and results with some known 
approximation schemes used for multivariate systems with colored noise. We shall show in particular that small 
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correlation time expansions cannot explain the anomalous diffusion exponent when truncated at any finite order. The 
partial summation of Fokker-Planck type terms, used to derive a 'best effective Fokker-Planck equation' for colored 
noise |171 Il8j leads to results that agree with ours. We emphasize that the noise considered in this work has a finite 
correlation time and its auto-correlation function does not have long time tails. 

This article is organized as follows. In section II, we analyse the case where the parametric fluctuations of the 
pendulum are modeled by Gaussian white noise, and explain heuristically why colored noise leads to anomalous 
scaling of the energy. In section III, we study the case of Ornstein-Uhlenbeck noise and explain how a recursive 
adiabatic elimination of the fast variable can be performed. This allows us to derive analytical results for the P.D.F. 
of the energy, which we validate with direct numerical simulations. In section IV, we compare our results with effective 
Fokker-Planck approaches. Concluding remarks are presented in section V. 



II. THE PENDULUM WITH PARAMETRIC NOISE 



The dynamics of a non-dissipative classical pendulum with parametric noise can be described by the following 
system of stochastic differential equations 

f2 = -(lu 2 + £(t)) sin0, (1) 

= O, (2) 

where 9 represents the angular displacement and Q the angular velocity. The energy E of the system is 

E= w 2 cos0. (3) 

2 v ; 

Equations (111211 describe the motion of a pendulum of frequency w whose suspension point is subject to a stochastic 
force proportional to the random function £(t). Our aim is to study how the stochastic properties of are transfered 
to the dynamical variables {9, Q) through a multiplicative and nonlinear coupling. We are chiefly interested in the 
case where £(t) has a non-zero correlation time, but we first consider Gaussian white noise. Eqs. Ill 121) . as well as all 
the stochastic differential equations below, are interpreted according to the rules of Stratonovich calculus. 



A. The white noise case 



The dynamics of an oscillator subject to multiplicative white noise has been studied by a number of investigators 
EH EU EH EU • In particular, Lindenberg et al. have shown J^] that the physical origin of the energetic instability 
of a stochastic oscillator can be quantitatively related to the parametric resonance of the underlying deterministic 
oscillator. In this section, an instability of the same type is found for the stochastic pendulum subject to parametric 
white noise, using the adiabatic averaging method. 

Let be a Gaussian white noise of zero mean value and amplitude T>: 

m = o, 

mm) = vs(t-f). (4) 

Using elementary dimensional analysis, we notice from Eq. that Cl *~ £, i.e., the angular velocity Q grows as 
t 1 / 2 and therefore 9 ~ i 3 / 2 . This observation can be put on a stronger basis by using the Fokker-Planck equation, 
associated with Eqs. JTJ and J5J), that describes the evolution of the Probability Distribution Function P t (9,Q): 

From Eq. J3J), we observe that the angular variable 9 varies rapidly as compared to Q. Thus, following |14| . we assume 
that, in the long time limit, the angle 9 is uniformly distributed over [0, 27r]. This allows us to average Eq. (0 over 
the angular variable and to derive an effective Fokker-Planck equation for the marginal distribution P t (Q): 

dP 1 = v& 2 P 1 

dt 4 an 2 ' u 



where we have replaced sin 2 9 by its mean value 1/2. From this effective Fokker-Planck equation, we readily deduce 
that f2 is a Gaussian variable with P.D.F. 



We thus recover that f2 grows as t x l 2 when t — > oo. Because E ~ f2 2 /2 (up to a term that remains bounded), we 
deduce the P.D.F. of the energy 



■nT>t V vt J 

From Eq. JHJ), we obtain the scaling behavior of the average energy 



=vJ=^W-l§Y (8) 



4 

and also the skewness and flatness factors 



(£> = ?*, (9) 



_ (£ 3 ) _ r(|) t$)W _ 5 ^ 

- ^5575- r(|)3/2 -^- 2.887..., (10) 
^) = m = ^ = ¥, (ID 



<£ 2 } 2 r(|) 

where T() is the Euler Gamma function. We conclude from Eq. © that the average energy of a frictionless pendulum 
with white parametric noise grows linearly with time. These results, Eqs. H9I11J) . agree with numerical simulations 
(see Fig.JTJ). 




FIG. 1: Stochastic pendulum with Gaussian white noise: Eqs. i l 1121 H are integrated numerically for V = 1. Ensemble averages 
are computed over 10 4 realizations. For numerical values of the pulsation uj = 1.0 and 0.0, we plot: (a) the average (E) and 
the ratio (E)/(Dt) (inset), (b) the skewness and flatness factors of E vs. time t. The asymptotic behavior of these observables 
agrees with Eqs. tSl Ull (dotted lines in the figures), irrespective of the value of u>. 



B. Scaling analysis for colored noise 

When the noise is colored, it is not possible to write a closed equation for the P.D.F. P t (9,il). However, the 
scaling behavior of the dynamical variables can be deduced from a self-consistent reasoning similar to that used in 
[THl ITfil | . Suppose a priori that we have the following scaling behavior in the long time limit 



(12) 



where a is an unknown exponent to be determined. We find from Eq. J5J that 9 ~ t a+1 , and Eq. Q can then be 
written as 



tl ~ £(*) sini 



Q + l 



(13) 



(We could have retained the deterministic term — up sin# but it would not affect this scaling analysis.) We now take 
£(i) to be a discrete dichotomous noise with correlation time r and with values ±1. The previous equation, then, 
becomes 



k=i 

We estimate the last integral by an integration by parts 



Q(t) - V e fe / sinx a+1 dx , with e k = ±1 . 

t_i J(k-l)r 
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t + T 
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Q + l 



-dx = 



cos(x Q+1 ) 



t+r 



+ 0(t' a - 1 ) 1 



and obtain 




sin a; 



Q + l 



dx 



(fe-l)r 



t/r 

E 



(fer) 2Q 



(14) 



(15) 



(16) 



This result is compatible with the a priori scaling Ansatz Ijl2|l only for a = 1/4. We thus conclude from this qualitative 
argument that, in the presence of colored noise, the energy E of the system defined in Eq. J2J grows as the square-root 
of time (as opposed to the linear growth obtained for white noise). A non-zero correlation time of the noise, however 
small, modifies the long time scaling behavior of the system. 

In the next section, we develop a method to put this qualitative analysis on a firm basis and derive precise analytic 
expressions that can be compared quantitatively with numerical results. 



III. THE AVERAGING METHOD FOR ORNSTEIN-UHLENBECK NOISE 



From now on, we consider the random noise £ to be an Ornstein-Uhlenbeck process, i.e., a Gaussian colored noise 
with correlation function given by: 



mm) 



^ e -l-'IA 
2r 



(17) 



where r is the correlation time of the noise. This noise £ can be generated from white noise via the Ornstein-Uhlenbeck 
equation 



1 



1 



(18) 



where rj(t) is a white noise of auto-correlation function T>S(t — t'). In the stationary limit, t,t' 3> t, the solution 
of Eq. (|18fl satisfies Eq. i|17|) . The pendulum with Ornstein-Uhlenbeck noise is thus written as a three-dimensional 
stochastic dynamical system coupled to a white noise rj(t): 



fl = -uj 2 sin 9 - £ sin 9 , 

6 = n, 



i 



i 



£ = --£ + -T»(t). 

r r 

The Fokker-Planck equation for the three-dimensional P.D.F. Pt(9, f2, £) is given by 



m 
at 
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((- 2 +Osin^) + 7|(^) + ^^ 



(19) 
(20) 

(21) 



(22) 



For colored noise, there is no closed Fokker-Planck equation for the original P.D.F. on phase space, Pt(9, fi), and only 
approximate evolution equations can be written. We shall not use any effective dynamics to derive our results but 
rather start with the exact three-dimensional Fokker-Planck equation Q22|l from which we shall integrate out the fast 
variable. 



A. Zeroth-order averaging 



We show here that the averaging procedure used in section Hi Al for white noise leads to erroneous results for colored 
noise. Averaging the Fokker-Planck equation ill'L'l) over the fast angular variable 9, we find that the marginal P.D.F. 
Pt(0, £) obeys the Ornstein-Uhlenbeck diffusion equation 

the variable f2 no more appears in this averaged Fokker-Planck equation and the associated stochastic two-dimensional 
system reads 

O = and £ = --£, + -v{t) . (24) 

T T 

In this averaged system, the angular velocity f2 is no more stochastic and, even worse, it is constant in time. By 
integrating out the fast variable 9 without taking into account the correlations between 9 and £, the dynamical variable 
Q, has been decoupled from the noise; in other words, the noise itself has been averaged out of the system. In fact, 
when the typical variation time of 9 {i.e. the time during which 9 varies by 27r) becomes less than r, the noise £ is 
roughly constant during a period of sm9. Thus, if 9 and £ are (wrongly) treated as independent variables, (sin#)£ 
is averaged to at the leading order. This problem did not occur in section III Al where £ was a white noise. In that 
case, the rapid fluctuations of the phase do not wipe out the noise, and (sin#) £ is averaged to C/V2 to yield Eq. JJJJ. 

In the next subsections, we develop an averaging scheme that allows us to eliminate adiabatically the fast variable 
while retaining the correlation terms. The idea is to define recursively a new set of dynamical variables that embodies 
the correlations order by order. This will enable us to derive sound asymptotic results for the pendulum with colored 
noise. In this scheme, Eq. I|24|) appears as a zeroth order approximation, and its correct interpretation is not that il 
is conserved but that its variations are slower than that of normal diffusion. 

B. First-order averaging 

Multiplying both sides of Eq. (|19|1 by f2, and using Eq. 1201) . we obtain 

fin = —Q (u 2 sin 9 + £ sin 6) = -u?9 sin 9 - £ (9 sin 9) . (25) 
Introducing the energy E of the system defined in Eq. J3J), this equation becomes 

f =i{* <*> 

We now transform the right hand side by writing it as a total derivative plus a correction term: 

^ = ^ cos9) - i cos9 . (27) 



Using Eq. 

we obtain 



d ,„ j. cos 6* . cos0 . „ N 

-(£:-£ cos0) = £ n. (28) 

at t t 



This leads us to define a new dynamical variable E\ 

Q2 

E x = E-£cos9= —~(lu 2 + f) cos 9 , (29) 

and to rewrite the stochastic system <|T§1 |2*T|) in terms of the set of variables (Ei,0,£): 

■ cos 9 . cos 9 , . ,„„ N 

Ei = £ v(t), (30) 



9 = n(E u 9,0 = ^2(Ei + {Loi+£,)cos9y (31) 
1 = --Z+-V(t). (32) 

T T 




FIG. 2: Stochastic pendulum with Ornstein-Uhlenbeck colored noise: Eqs. 119H2H are integrated numerically for u) = 1.0 and 
(V, t) = (1, 0.1), (1, 1), (1, 10), (10, 0.1), (10, 1), (10, 10), (100, 0.1), (100, 1), (100, 10). Ensemble averages are computed over 10 4 
realizations. We plot the ratio (E 2 )/(T>t/T 2 ) vs. time t. The dotted line in the figure corresponds to the asymptotic value 0.25 
predicted by second-order averaging (Eq. 155^ 1. Note that convergence is slower for smaller values of the correlation time r. 



The advantage of this system as compared to the previous one H19ll20ll2"TJ) is that the white noise rj(t) now appears 
in the equation for the dynamical variable E\ and this white noise contribution will survive the averaging process. 
The Fokker-Planck equation for the P.D.F. P t (E 1 ,6,^) associated with Eqs. $M EH and 03 reads: 

Averaging this equation with respect to 8 leads to the following evolution equation for the marginal distribution 

dt ~ r OH V n ) + 4r 2 8El + 2r 2 ' 1 ' 

The variable E\ now appears in the averaged Fokker-Planck equation, and the associated stochastic two-dimensional 
system reads 

Ei = ^i(i), (35) 
i = --£+-»?(t), (36) 

T T 

where rj(t) and rj\(t) are two independent white noises of amplitude T>. This effective dynamics is exactly solvable 
(the variables E% and ^ are decoupled thanks to the absence of the second-order cross-derivative term in Eq. I|34|l). 
If we compare this system with the one obtained by naive averaging l)24|). we observe that the dynamical variable 
Ei (and therefore the energy of the system) is no longer constant, rather it grows as the square-root of time. When 
t — * oo, Ei becomes identical to the energy of the system (up to terms that remain finite). We therefore determine 
the long time statistics of the energy from Eq. (|35|l by identifying E\ to E and by imposing the physical condition 
E > 0. The energy is thus a Wiener process on a half-line and its P.D.F. is given by 

P t (£). jL=ex P (-^) with E>0. (37) 

From this expression we calculate the first two moments of the energy 

w _£(*)' aMM (*)\ m 

<E 2 > = \™- (39) 



These expressions provide scaling relations between the averages, the time t and the dimensional parameters of the 
problem, T> and r. We have verified numerically that these scalings are correct (see Fig. However, the prefactors 



that appear in Eqs. 1(38(1 and 1(39(1 are pure numbers and do not agree with the results of our numerical simulations. 
We conclude that Eq. ((38(1 is exact at leading order as it gives the correct asymptotic scaling for the energy, E oc i 1 / 2 
but fails to provide the prefactors. The reason is that some correlations between 9 and £ have still been neglected in 
the averaging procedure. Carrying out the calculations to the next higher order will enable us to derive the correct 
expressions for the prefactors. 



C. Second-order averaging 



More precise results can indeed be derived by applying recursively the procedure described above. In the Langevin 
equation 1(30(1 for E\, we perform one more 'integration by parts' and obtain 

• fcosO 9 cos 9 , . d /£sin6>\ sin 9 • f2 sin# . cos 9 , . , ,„ N 

* = V n - — " {t) = dt ( V ) - ^ + w — « " — • (40) 

This leads us to introduce a new variable E2 defined as 



£sin6» n 2 , 9 ^ „ fsin9 
^ = ^-'— = --^+0^9-'—. (41) 



Using Eqs. JTSJ and (|3TJ, Eq. igDJ becomes 

• £sin6> £sm 2 9 , 9 _ / „ sin#\ rcfi) , 

£2 = " W - + " ( cos 9 + ^n)r- ^ 

In this equation we must express the variable fl in terms of E2, 9 and £. Inverting the relation (|41|) . we deduce that 

n = (2^)4 + (. 2 + O^r + ^1 + o f (43) 

(2£ 2 )2 r(2ii 2 ) V-^j / 
where we have retained terms up to the order 1/E 2 . From Eq. ((42(1 . we deduce the Langevin equation for £?2 

£2 = J E (E 2> e,0+D E (E 2 ,6 > Z) 1 ^ +0 (4r) , (44) 



where we have defined 



D E (E 2 ,e,£) = ~(cos9+ ^-r ) . (46) 



r(2 J E 2 ). 

This equation, combined with Eqs. I|20l and 121(1 . defines a three-dimensional stochastic system for the variables 
(E 2 ,6,£). The Fokker-Planck equation for the P.D.F. P t (E 2 ,9,£) is 

F)P r) r) 1 r) 

T) ( 8 r) f) 2 cl BP r) 2 P 1 

+ 2^ {^W^ + + 8E 2 D ^ + (4?) 

We now integrate out the fast angular variable 9 from Eq. 1(47(1 . retaining only the leading term in the average of the 
expression -^-D E j^-{D E Pt), (recalling that d/dE 2 scales as E^ 1 , the contribution of the subdominant terms is of 

the order of E 2 and is negligible in the long time limit). We thus obtain the following evolution equation for the 
marginal distribution P t (E 2 ,£): 

dP t d fu 2 £, + £, 2 ~\ 1 d V (id 2 P t , d 2 P t ^ 

-Pt 



dt dE 2 \ AtE 2 7 rd£ V s V ' 2r 2 2 dE 2 di 



Although the cross-derivative terms between E2 and £ vanish, the two variables are coupled through the drift term. 
From Eq. H48J1 we derive an effective, two-dimensional, stochastic system in E2 and £: 



E > = -^T + ^ 2(<) ' (49) 

£ = --Z+-r){t), (50) 
r r 

•q{t) and T72 (t) being two independent white noises of amplitude T>. If we compare this system deduced by a second- 
order averaging with the one obtained at first order Eqs. (|35lrt6[) . we observe that the equation for the dynamical 
variable E2 contains, besides a noise term, an effective potential that scales like I/E2 and that involves the other 
variable £. This effective potential diverges at Eq, = and constrains the energy to stay positive: we do not need 
anymore to impose this physical condition arbitrarily. 



D. Analytical results 

In the coupled system I|49I50I) , the fast angular variable has been eliminated and the dimensionality of the original 
problem has been reduced by one. However, Eq. 1|49[) is nonlinear and is not exactly solvable. Nevertheless, the long 
time behavior of E2 can be deduced from the following reasoning. Recalling that £ 2 has a finite mean value, equal to 
T>/2t, we rewrite the evolution equation (|49|l of Eq as 

+ ^» - 55 <■-«*»■ < 51 > 

This Langevin equation contains two independent noise contributions : a white noise 772 (t) , and a (non Gaussian) 
colored noise, (w 2 £ + £ 2 — (£ 2 )), of zero mean value and of finite variance. This colored noise is multiplied by a 
prefactor proportional to I/E2 and, because E2 goes to infinity with time, it becomes negligible in the large time 
limit in comparison with the white noise term. Thus, Eq. I|51f) reduces asymptotically to 

h = -jE 2 + Wr m(t) - (52) 

The variables E2 and £ are decoupled at large times: the effective problem is thus one-dimensional. In the long time 
limit, the variable E2 is identical to the energy E up to finite terms. We thus obtain the asymptotic P.D.F. of the 
pendulum's energy by explicitely solving the Fokker-Planck associated with Eq. il~>L't 

PJE) = —^L—E-i exp (-—) . (53) 

v ; r(i)(pt) 1 /4 y \ vt ) K ' 

Hence, E is not a Wiener process on a half line: its P.D.F. is not a simple Gaussian. From Eq. \b'6\. we calculate the 
first two moments of the energy 

V2tt [Vt\ 1/2 [Vt^ 1/2 



n\) 2 



( E ) = ^7T^ I ~2 ) * °- 338 ( -2 ) > ( 54 ) 



^) - ~. (55) 
Besides the skewness and the flatness factors are 

SW-ESpi.^* 2.028..., (56) 

F(E) = = 5. (57) 

The functional dependence of the moments on time t and on the parameters T> and r is the same as that obtained in 
section UlI Bl But the prefactors, which are absolute numbers, are different (compare Eqs. (j38l39|l with Eqs. 



Numerical simulations of the dynamical equations (jl9l 1201 and !21|l . shown in Figs. Eland |3| agree quantitatively with 




FIG. 3: Stochastic pendulum with Ornstein-Uhlenbeck colored noise: Eqs. 119H2H are integrated numerically for T> = 1, 
t = 1. Ensemble averages are computed over 10 4 realizations. For numerical values of the pulsation lo = 1.0 and 0.0, we plot: 
(a) the average (E) and the ratio (E) /(Dt/r 2 ) 1 / 2 (inset), (b) the skewness and flatness factors of E vs. time t. The asymptotic 
behavior of all observables presented agrees with Eqs. 1)54145 7|l (dotted lines in the figures), irrespective of the value of to. 

the predictions of Eqs. (|54I57(1 . In particular, we notice that the asymptotic P.D.F. P t (E), and therefore the moments 
of the energy, do not depend on the value of the mean frequency to (see Fig. |2J). 

Our averaging technique thus provides sound asymptotic results for the energy of the stochastic pendulum: this 
technique not only yields the correct scalings but also leads to analytical formulae for the large time behavior of the 
energy. This averaging method could be carried over to the third order to calculate the subdominant corrections to 
the P.D.F. of the energy. However, we shall not pursue this course any further: the calculations become very unwieldy 
and the agreement between the analytical results and the numerical computations is already very satisfactory at the 
second order. 

We emphasize that the white noise and the colored noise cases fall in two distinct universality classes because the 
long time scaling exponents are different. The colored noise scaling will always be observed after a sufficiently long 
time provided that the correlation time r is non-zero. However, the effect of this correlation time appears only when 
the period T of the pendulum is less than r. This period, which is proportional to fi , decreases with time. At 
short times, the period T is much greater than r and white noise scalings are observed. At large times, T < t and 
colored noise scalings are satisfied. The crossover time t c is reached when T ~ ~ r; using Eq. 10, which is valid 
for t < t c , we obtain 

t c ~ (Vt 2 )- 1 . (58) 

Hence, when the correlation time r becomes vanishingly small, the crossover time diverges to infinity and the colored 
noise regime is not reached (simulation times needed to observe the colored noise scalings become increasingly long). 
In Fig. 21 we plot the behavior of the mean energy vs time for T> — 1.0 and r = 0.01, 0.1 and 1. For r = 1, the colored 
noise scaling regime is obtained from the very beginning and the curve for (E) has a slope 1/2 in log-log scale. For 
t = 0.1, at short times (E) oc t whereas at long times (E) oc t 1 ^ 2 . The crossover is observed around t c ~ 100 = t~ 2 . 
For r = 0.01, the curve for (E) has a slope 1 and the colored noise regime is not reached in this simulation (t c ~ 10 4 ). 

The averaging method provides analytical results for t^> t c and f<t c . The intermediate time regime is described 
by a crossover function of the scaled variable t/t c |15| . which cannot be analyzed by our technique. In the next 
section, we compare our method with two well-known approximations based on effective colored noise Fokker-Planck 
equations. One of these approaches will enable us to derive approximate formulae for t c and for the crossover function. 

IV. COMPARISON WITH OTHER APPROXIMATION SCHEMES 

The main difficulty for the study of a Langevin equation with colored noise stems from its non-Markovian character 
PJ: there exists no closed Fokker-Planck equation that describes the evolution of the P.D.F. of the dynamical variables. 
The process can be embedded into a Markov system if the colored noise itself is treated as a variable. However, this 
mathematical trick increases the dimension of the problem by one and the noise has to be integrated out in the end: 
exact calculations can be carried out with this method only in the case of linear problems. 
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FIG. 4: Stochastic pendulum with Ornstein-Uhlenbeck colored noise: Eqs. I)19H21|I are integrated numerically for ui = 0.0, 
V — 1.0 and r = 0.01,0.1 and 1. Ensemble averages are computed over 10 4 realizations. We plot the average energy vs. 
time. The crossover between white noise behavior ((E) ~ t) and colored noise behavior ((E) ~ t 1 ^ 2 ) occurs later for smaller 
correlation time r. 

A colored noise master equation can be rigorously derived (e.g., with the help of functional methods) but it involves 
correlations between the dynamical variables and the noise "231 . The equation of motion for these correlations 
involves higher order correlations, and so on. Since this hierarchy must be stopped at some stage, this question is a 
genuine closure problem. Many different app roaches have been devoted to derive effective Fokker-Planck equations, 
such as short correlation time expansions Il7tli~8l . l25l . l26l ]. unified colored noise approximation ~~~|, projection methods 
[28| and self-consistent decoupling Ansatze 29] (for a general review see 30]). 

In section IIIII to derive the long time behavior of the stochastic pendulum with colored noise, we did not use any 
closure approximation but started from the exact Fokker-Planck equation for the system, treating the noise as an 
auxiliary variable. Thus, we did not make any hypothesis on r and our analytical formulae are valid for any value of 
the correlation time (for t larger than the crossover time given in Eq. """J). In this section, we compare our results 
with those that can be derived from some well-known approximation schemes. 

A. Small correlation time expansion 

An approximate evolution equation for the P.D.F. of a Langevin equation with colored noise can be derived in 
the case of short correlation times (i.e., in the white noise limit) by expanding the colored noise master equation 
around the Markovian point |26j, |3Tj • This procedure leads to a Fokker-Planck type equation, with effective drift and 
diffusion coefficients. Applying to our system (|19M20ll2~"|l the small r expansion derived in [2f| for arbitrary stochastic 
equations with colored noise, we obtain, at first order in r, the effective Fokker-Planck equation for Pt(0, SY) 

8P t 32?r dP t , £>(1-Tft) ■ 2a d 2 P t , Vt d 2 ,. 2flp x 

-dT = -de {m) -— smdcosd m + 2 sm ^ + -^ sm flp *)" (59) 

To simplify the discussion, we have taken the mean frequency u> equal to zero. After integrating out the rapid 
variations of 6 an averaged Fokker-Planck equation is obtained which is similar to Eq. "dJ), but with an effective 
diffusion given by T> c s = X>(1 — rfi) (because T> e g becomes negative for large values of SI, Eq. """} is valid only over 
the restricted region of positive diffusion). From dimensional analysis, we observe that such an effective Fokker-Planck 
equation leads to a normal diffusive behavior, SI ~ i 1 / 2 and E ~ t, and therefore cannot account for the results we 
have obtained. 



B. Best Effective Fokker-Planck equation 

We now turn to another small r approximation, which intends to improve the first-order effective Fokker-Planck 
equation Ij59(l by summing contributions of the type T>r n (where n is an integer larger than 1). The resulting equation 
has been christened 'Best Fokker-Planck Equation' (B.F.P.E.) by its proponents Il7lll8ll3l| . Although this approach 
is not free from drawbacks and is known to lead in some cases to unphysical results |32|. we show that for the system 
studied in this work, the B.F.P.E. leads to results that agree with ours. 



An approximate evolution equation for the P.D.F. Pf(Q,6) is given by the second-order cumulant expansion |l| of 
the (stochastic) Liouville equation associated with Eqs. I|19l and l2~n|l : 



dP t 



L Q P t + [ dx(L 1 (t)aq){L x)L 1 (t-x)exp(-L x))Pt (60) 
Jo 

where the differential operators are defined as 

LoPt = ~§8 imh (61) 

U{t)Pt = -^(C(t)ainflP t ). (62) 

(Here again we take lo = 0). In the Appendix, we evaluate the right hand side of Eq. (|6Uf) and derive, in the limit 
t — ► oo, the following B.F.P.E. for the classical pendulum with colored multiplicative noise 

dP t d . V d 2 ( sin 2 6 ~ rtt sin 6 cos 6 \ Vrd „ d ( (1 - (tVi) 2 ) sin 9 - 2rfl cos 6 \ , s 

1* = + 2 W { i + ijW Pt ) + -on smd de { ' ) ' ( } 

Integrating out the fast variable 6, we obtain an averaged B.F.P.E for Pt(O), the probability distribution of the slow 
variable, 



dP t _ V d / 1 dP t \ 

~~dt ~ 1WL\i + {T9) 2 ~m) 



(64) 



For t = 0, this equation is identical to the averaged white noise Fokker-Planck equation JHJ). For a non-zero correlation 
time, this equation predicts correctly that grows as t 1 / 4 : this is straightforward from scaling. The crossover between 
white and colored noises is observed for £1 ~ 1/r and this is consistent with the discussion that led to Eq. I|58|l . 
Although Pt(O) is not explicitely calculable, Eq. (|64Jl implies the following identity 



^V) + i(tt 2 H^. (65) 



When t is small, the quadratic term dominates over the quartic term, and we recover (SI 2 ) = 2(E) ~ yt, in agreement 
with Eq.Q. When t is large, the quartic term is dominant and one deduces that 



(E 2 ) = -(n 4 )^--t. (66) 



4 X 1 At 2 

This result is identical to our Eq. (|55|l . which was validated by numerical results (see Fig.|2Jl. The identity Ij65(l can 
also be used to derive an approximate scaling function for the mean energy. Let us define the flatness (f> of Pt (ft) as 

(r> 4 ) = (j> (n 2 ) 2 . (67) 

Rigorously speaking is a function of time, but it remains a number of order 1. For simplicity, let us assume that <p 
is constant. Substituting Eq. I|67|l in Eq. (|65|l . and solving for (ft 2 ), we obtain 



This scaling function explicitely describes the evolution of the energy of the oscillator as a function of time. It contains, 
in particular, the linear behavior at short time and the t 1 / 2 growth at large time. The crossover between these two 
scaling regimes occurs when Dr 2 t ~ 1, i.e., precisely at the crossover time given in Eq. (|58() . 

The B.F.P.E. approach has thus allowed us to derive short and long time scalings and to understand qualitatively 
the evolution of the system at intermediate times. The B.F.P.E. is derived from a second-order cumulant expansion 
in which higher order terms, that are not Fokker-Planck-like, are discarded. This approximation has been criticized 
|32| because the neglected higher-derivative terms can be of the same order as the terms that have been retained in 
the B.F.P.E. However, for the stochastic pendulum, these neglected terms do not change qualitatively the large time 
behavior of the dynamical variables. In contrast, the averaging technique that we have generalized here is not based on 
any a priori expansion and provides reliable results at least for the stochastic system studied in this work. Of course, 
the approach advocated here has been developed for one particular problem and does not have the generality and 
versatility of effective Fokker-Planck methods. Nevertheless, we strongly believe that the recursive averaging scheme 
developed here for the stochastic pendulum can be extended to other nonlinear one-dimensional systems subject to 
multiplicative or additive colored noise |33| . 



V. CONCLUSION 



A nonlinear pendulum subject to parametric noise undergoes a noise-induced diffusion in phase space. The charac- 
teristics of this motion depend on the nature of the randomness: when the noise is white the energy of the pendulum 
grows linearly with time, whereas it varies only as the square root of time when the noise is colored. This change of 
behavior is due to destructive interference between the displacement of the pendulum and the noise term: the effect of 
the colored noise is partially averaged out by fast angular variations as soon as the period of the pendulum becomes 
smaller than the correlation time r of the noise. We have carried out an analytical study of this model by defining 
recursively new coordinates in the phase space and averaging out the fast angular variable. At zeroth order, the only 
information obtained is that the motion is subdiffusivc; at first order, this procedure provides the correct scalings; at 
second order, a quantitative agreement with numerical simulations is reached. 

We emphasize that our method is different from the usual approximations that involve effective Fokker-Planck 
equations in which the colored noise does not appear as an auxiliary variable. Our averaging procedure integrates out 
the fast dynamical variable and leads to an effective stochastic dynamics for the slow variables with colored noise. 
Whereas usual effective Fokker-Planck equations are valid only for small noise and for short correlation times, we 
do not make any hypothesis on the amplitude of the noise or on the correlation time r. However, the asymptotic 
subdiffusivc regime is reached earlier for larger values of t. Our results agree with those derived after averaging 
the 'Best Fokker-Planck Equation' (B.F.P.E.) approximation (which is obtained from a summation of a truncated 
cumulant expansion of the Liouville equation) . This approximation is also useful to draw a qualitative physical picture 
of the system and allows to calculate approximate crossover functions between short time (white noise) and long time 
(colored noise) regimes. 

The recursive averaging scheme that we have used here for the stochastic pendulum can be extended to other 
systems coupled with colored noise. In particular, we believe that thanks to this method a precise mathematical 
analysis of the long time behavior of a nonlinear oscillator subject to multiplicative or additive colored noise can be 
carried out |3^| . 
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APPENDIX A: DERIVATION OF EQ. (l63t 



In this appendix, we derive the B.F.P.E. (Eq. I|63|) ) following the procedure of 0,0]- We evaluate the right hand 
side of Eq. (jfit)|l by applying the following operator formula 

exp(A)i?exp(-A) = B + [A, B] + |L4, [A, B}] + ±[A, [A, [A, B}]} + . . . , (Al) 

with A = Jjq and B = ~Li(t — x). We thus have to calculate some commutators of the two operators Lo and Li(i — x) 
defined in Eqs. I|tjl|) and Hti2[l . By induction, we derive the following expression for the n-th commutator 

T n = [L , [. . . , [L 0) [LoMit x)}} ...}=Z(t- x) (J^H[ n \n,6) + ^^(fi,*)) , (A2) 
where the functions and satisfy the recursion relations 

o rr{n — 1) 

B§> = H?- 1] ~ . (A4) 



The first few terms can be calculated explicitely and we obtain 





= — sin 8 


and 


iff 


= 0, 




= il cos 8 


and 


H p 


= — sin 8 


iff 


= fl 2 sin 8 


and 




= 2fl cos 



(A5) 
(A6) 
(A7) 



The general solution for the recursion ljA4fl is readily found: 

H{ n) = (-lf-^sin^ + n-), (A8) 



H [ 2 n> = nH^ 1 ' = (-l) n - 1 nn n ' 1 cos{9 + n-). (A9) 



From Eqs. (|A1|) and i|A2(l . we deduce the following identity 



x" 



(Li(t)exp(L x)Li(t - cc)exp(-L x)) = J2 n =o ~ r( L i(*) T ' 



= ~£i'**0 ££o ^(m(t-x))(^Hi n \a t 0) + ^Hi n \^ . (A10) 

Substituting in this equation the expressions l)A9|l of Hi and H2 ; an d the autocorrelation function (|17|l of the 
Ornstein-Uhlenbeck noise, we find that the right hand side of Eq. is given by 

jf dx{Li{t)exp(L x)Li{t-x)exp(-L oX ))=- — — sm6 \^Hi{£l,0,t) + —Hi{£l,0,t)\ , (All) 
where we have defined 

HiMt) = g i" ) = E ° „, (-lr-V^O^sin^ + n^), (A12) 

ra=0 ' n=0 

H 2 (n, 9,t) = Y joj^£^I g W = y f^^-' ^^-i^y.-! cos(0 + n 7r } _ (A13) 

n=0 n=0 

In the limit i — > 00, the integral f *^ T dxx n e~ x converges to nl, and the series defining Hi and H2 can be calculated 
in a closed form. We finally obtain 

„ , ,^ n ^ sin# — ril cos 9 . , „ ,^ 

Hi (fl, 0,oo) - -r 1 - ^ , (A14) 

„ . 9 (1- (tIT) 2 ) sin (9- 2rO cos (9 . 

n 2 (n,e,oo) = t - — ^ — ^ 5 . A15 

(l + (rft) 2 ) 2 

This completes the proof of Eq. (|rJ5|) . 
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